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Abstract. We perform three-dimensional hydrodynamical sim- 
ulations of the coalescence of binary neutron stars. We include 
the emission and backreaction of gravitational waves into the 
Newtonian "Piecewise Parabolic Method". The use of the phys- 
ical equation of state (EOS) of Lattimer & Swesty (1991) allows 
us to take into account the production of neutrinos. We evaluate 
our models for the efficiency of vv annihilation in the surround- 
ings of the coalescing neutron stars. The corresponding energy 
deposition prior to and during merging turns out to be 2-3 or- 
ders of magnitude too small to power a typical 7-ray burst with 
an energy output of ~ (10 51 /47r) erg/sterad at cosmological 
distances. Analytical estimates of the subsequent evolution of 
the disk which possibly surrounds the central black hole show 
that even under the most favorable conditions the energy pro- 
vided by vv -> e~e + -> 77 falls short by at least an order of 
magnitude. 
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1. Introduction 

The distribution of 1 1 22 observed 7-ray bursts as reported by the 
BATSE 3B Catalog is isotropic over the sky and non-uniform 
in distance (Meegan et al. 1995). One of the currently most 
popular models for the origin of the 7 emission is the merger of 
two neutron stars (NS-NS) or of a neutron star and a black hole 
(NS-BH) at cosmological distances (Paczynski 1986, Goodman 
1986,Eichleretal. 1989,Narayanetal. 1991, Paczynski 1991). 
These would occur at a rate of 10~ 4 -10~ 5 yr -1 per galaxy and 
thus with a sufficiently large frequency to explain the observed 
burst rate (Phinney 1991, Narayan et al. 1991). Only about 1% 
of the several 10 53 erg of energy released in the merging would 
be enough to power a burst, and the compactness of the source 
might account for the variations in the 7-ray flux on timescales 
as short as 1 ms. 

The non-thermal 7-rays are supposed to be produced by a 
fireball of relativistic e~e + pairs that are created by the anni- 
hilation of vv pairs in the vicinity of the hot, merged object 
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(Goodman et al. 1987). However, the large v fluxes also drive 
a baryonic wind by energy deposition near the neutrinosphere 
(Duncan et al. 1986, Woosley & Baron 1992) which "pollutes" 
the cloud of 7's and e~e + pairs. The baryonic load prevents 
relativistic expansion and keeps the flow optically thick, lead- 
ing to a degradation of the 7-rays. Thus the energy release will 
be mainly kinetic rather than in photons (Paczynski 1990). If 
spherical symmetry is broken, as in the case of the merger sce- 
nario where a rapidly rotating disk can be formed and material 
is pulled away from the rotation axis by centrifugal forces, then 
a baryon-free funnel along the rotation axis may allow relativis- 
tic beams of 7's and e~e + to escape (Meszaros & Rees 1992a, 
Mochkovitch et al. 1993, Woosley 1993). Annihilation of vv 
along the system axis might lead to matter ejection with Lorentz 
factors of T ~ 10 2 -10 3 (Mochkovitch et al. 1995) which are in 
the right range to enable copious 7 production during shock in- 
teraction with ambient interstellar gas (Rees & Meszaros 1992). 

In this Letter we present 3D hydrodynamical simulations of 
the merging of binary NSs using the physical EOS of Lattimer 
& Swesty (1991) and including the effects of the backreaction 
of gravitational waves on the hydrodynamic flow (see Ruffert 
et al. 1996). The loss of lepton number and energy by j/'s is 
simulated with an elaborate leakage scheme that was calibrated 
using diffusion results. It takes into account the emission of all 
types of z/'s by thermal processes and of v e and v e also by the 
/3-processes. Neutrino opacities and diffusion timescales are 
calculated from v-n,p scattering and the inverse /3-reactions 
(Ruffert et al. 1996). We evaluate our models for the efficiency 
of vv annihilation and investigate the question whether the 
corresponding energy deposition is powerful enough to account 
for the energetics of burst events at cosmological distances. 
That is the crucial question when associating 7-ray bursts with 
NS-NS or NS-BH mergers, independent of the detailed, very 
complex processes which lead to the emission of observed 7- 
rays. 

2. Thermal evolution and neutrino emission 

Pre-merging evolution. Kochanek (1992) showed that un- 
less the volume-averaged shear viscosity is larger than r\ ^ 
1 29 g cm" 1 s ~ 1 the structure and evolution of the NS binary sys- 
tem is not significantly affected by viscous effects. Tidal locking 
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(synchronization) of NS-NS orNS-BH systems just before tidal 
disruption requires a dynamic viscosity r) J> 10 31 gcm _1 s _1 
(Kochanek 1992, Bildsten & Cutler 1992). 

The microscopic shear viscosity of NS matter is much 
smaller. For matter hotter than the superfluid-transition tem- 
perature (T <; 10 9 K) it is dominated by neutron-neutron in- 
teractions. For the density range 10 14 & p Ss 4 • 10 14 g/cm 3 it 
was calculated by Flowers & Itoh (1979) and can be fitted to an 
accuracy of a few percent as r] n = 1.1 • \0 16 p\^ /T$ g cm _1 s _1 
where p l4 = p/10 14 g/cm 3 and T 9 = T/10 9 K. For NSs cooler 
than 10 9 K the neutrons become superfluid and their effective 
scattering cross section with other particles quickly falls to zero. 
The dynamic shear viscosity is then dominated by electron- 
electron scattering and is approximated at NS densities to within 
a few percent by rj e = 6 • 10 16 (pi4/T9) 2 gcm _1 s _1 (Cutler & 
Lindblom 1987). At temperatures of a few 10 9 K NSs become 
opaque for v's. Thus shear viscosity as a result of momentum 
transport due to the scattering of v's is potentially important. 
Assuming as main opacity source the scattering of v off n,p 
(and the charged current absorptions for v e and z> e ), and tak- 
ing into account the contributions of 3 types of nondegenerate 
(zero chemical potential) z/s and i>'s in local thermodynamic 
equilibrium with the stellar matter, the result of van den Horn 
& van Weert (1981) for the v shear viscosity can be generalized 
toi7„ ~ 8.6 • lO^/pugcm-^- 1 . 

Reisenegger & Goldreich (1994) showed that the energy 
absorbed by the coalescing NSs is only ~ 10~ 6 of the orbital 
energy even for stellar oscillation modes that are resonantly ex- 
cited by the strong tidal forces when the two stars come close. 
This was confirmed by Lai (1994) who also found that tidal 
dissipation heats up the stars primarily by resonant excitations 
of the g-modes. Using the microscopic viscosity he found that 
the core temperature is increased from less than ~ 10 6 K in NSs 
older than 10 8 yr to only ~ 10 8 K before the stars come into 
contact. To obtain a temperature of about 10 11 K as estimated 
by Meszaros & Rees (1992b) would require a viscosity that is 
larger than the microscopic one by a factor of 10 10 -10 12 , cor- 
responding to a dynamic viscosity near the maximum possible 
value in the NS, r? max « pcR ns « 10 31 gcm _1 s _1 (R DS « 10km 
is the NS radius). Such a high value of the viscosity appears 
rather unlikely and the NSs should stay comparatively cool 
before they merge. 

Merging - numerical simulations. We performed hydrodynam- 
ical simulations of the NS merging for three different setups of 
the initial NS spins. In one case both NS spins were parallel 
to the orbital angular momentum vector (model B64), meaning 
rigid rotation of the NS-NS system, in another case both spins 
were in the opposite direction (C64), and in a third case the NSs 
had no spins (A64). The two identical, cool Newtonian NSs had 
a baryonic mass of M ns « 1 .63 M© each (corresponding to a 
gravitational mass of 1.4-1.5 M©, dependent on the EOS), a 
radius of i? ns « 15 km, and spin periods being the same as the 
initial orbital period. Initially the stars are spherical and orbit 
around the common center of mass with a center-to-center dis- 



tance of 42 km and an orbital period of 2.6 ms. The orbit decays 
due to gravitational wave emission, dynamical instability sets 
in, and within ~ 1 .5 ms the NSs start to merge. 

We chose an initial, central temperature in the stars around 
kT ~ 7 MeV, corresponding to a thermal energy of roughly 
3% of the degeneracy energy for given density p and electron 
fraction Y e . As the stars approach each other, they expand due 
to tidal stretching. This leads to a slight drop of the core tem- 
peratures to about 6 MeV at the centers. These values are close 
to 10 11 K and are therefore of the order of the estimates of 
Meszaros & Rees (1992b) for the final stages of the binary be- 
fore coalescence. In view of the discussions above we consider 
them to be far too large. Even if the extreme shear dissipa- 
tion of a Keplerian disk were present in the coalescing NSs 
(certainly a strong overestimation of the true shear motions) 
and no v cooling took place, the heating to a temperature of 
~ 6 MeV during the last ~ 24 s of the binary spiral-in (start- 
ing with cold stars at a separation of 20i? ns ) would require a 
constant viscosity 77 ~ 10 26 gcm _1 s _1 which is larger than all 
microscopic viscosities by at least 4 orders of magnitude. The 
assumed high initial temperatures, however, favor copious and 
early v emission and thus serve the directions of investigation 
in this work. 

During merging the stars heat up drastically because of 
dynamical compression and dissipation of kinetic energy. In 
particular at the shear interface where the two stars come into 
contact first, temperatures between 1 and 30 MeV are obtained, 
dependent on the NS spins. After ~ 3ms from the start of the 
simulations, the two NSs have fused to a compact object with 
mass M ss 3 M Q , a radius of ~ 20 km, and densities p > 
(1 ... 2) • 10 14 g/cm 3 . The central body reaches temperatures of 
more than 50 MeV. It is surrounded by a thick disk with orbital 
velocities of 0.3-0. 5c, densities p ~ 10 12 g/cm 3 , and a mass 
Miisk ~ 0. 1-0.2 Mq. During the first few ms after merging this 
disk is heated from initially 6 MeV to ~ 10 MeV by dissipation 
and compression in strong pressure waves that are created by the 
wobbling and oscillations of the central object. Since the EOS 
used cannot stabilize NSs more massive than ~ 2.2 M©, we 
expect the massive object to collapse to a BH within a few ms. 
About 10% of the initial angular momentum are carried away 
by gravitational waves. With a relativistic rotation parameter 
a = Jc/(GM) <, 0.5-0.7 which is smaller than the critical 
value, a cr it « 1 (Stark & Piran 1985), the central body cannot 
be stabilized by rotation. Only matter that has a sufficiently 
large specific angular momentum, RiV^ ^ rv ^ R\c, where 
Ri ss 6GM/c 2 is the radius of the innermost stable circular orbit 
and i> k = ^jGMfRi the Newtonian Keplerian velocity at radius 
r = Ri, might remain in a disk around the BH. We estimate that 
with a typical rotation velocity v « 0.25c only gas with orbits 
r ^ 44 ... 1 07 km will have a chance not to be swallowed up by 
the BH on a dynamical timescale. In our calculations we find 
that at most 0.01-0. 15 M© fulfill this condition, most likely less. 
This value is extremely sensitive to the angular momentum in 
the NS-NS system and thus to the assumed NS spins. About 
10-30% of the initial angular momentum might be stored in the 
disk. 
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Fig. 2. Average energies of emitted i/ e , v e , and v x for the three models 



The luminosities of v e , v e , and v x (the latter representing 
the sum of the contributions of v^, i> M , v T , and v T ), are plotted 
in Fig. 1 for all 3 models. The maximum total v flux of 1- 
1 .5 • 1 53 erg /s is reached later than ~ 6 ms and is about 2 orders 
of magnitude smaller than the peak luminosity in gravitational 
waves which is obtained during the dynamical merging of the 
two stars (at ~ 1 .5 ms). The total energy emitted in gravitational 
waves is 0.01-0.02 M Q c 2 , while z/s radiate away an energy of 
only 2-3 • 10~ 4 M Q c 2 within about 10 ms after the start of the 
simulations. The corresponding average energies of the emitted 
z/s are shown in Fig. 2. The emission is dominated by v e with 
a mean energy of 18-22 MeV. They have a flux of « (5.5- 
6.5) • 10 52 erg/s in the stationary situation for t ^ 6 ms. The 
luminosities of v e and v x are smaller, L Ve : : 4L„ m ss 
1:4:3, and their mean energies are (e„ e ) ss 12-13 MeV and 



{e Vx ) ss 26-28 MeV, respectively. By far most of the v emission 
comes from the disk. After having reached its maximum interior 
temperatures kT « 10 MeV at t *S 5 ms, the disk emits more 
than 90% of the total flux. While the very dense, much more 
opaque core loses energy by v emission at a typical rate of 
~ 10 31 -10 32 ergcm 3 s \ the disk emits z/s at a rate of (3- 
10) • 10 32 erg cm _3 s _1 . When decompressed and heated, n-rich 
(Ye ~ 0.01-0.09) NS matter becomes moderately degenerate 
and the /3-equilibrium is shifted to fi e + fi p — fi n = //„ e < 0. 
In such a situation the process e + + n — > p + v e is faster than 
e~ + p -> n + v e and leads to the dominant production of v e in 
the disk. 
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fluxes of v e , v e , and v x = 4i/ u for the three models 7- 
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Fig. 3. Integrated vv energy deposition rate in the region around the 
merged neutron stars out to an equatorial radius r in the vertical range 
< \z\ < 64km and in the range 20km < |z| < 64 km, respectively 



Our models were evaluated for the local energy deposition 
rates by vv annihilation in the vicinity of the merged NSs in 
regions where p < 10 11 g/cm 3 and where the energy loss rate 
by v emission is less than 10 30 erg s _1 cm -3 (for details see 
Ruffert et al., in preparation). Fig. 3 displays the cumulative 
energy deposition rate dE vD (r)/dt, integrated out to a radius 
r in the equatorial plane (r < 64 km), once for the vertical 
range \z\ < 64km, and another time only for the volume that 
lies above and below the thick disk around the central object 
(20km < \z\ < 64km). For the disk geometry of our models 
we find with the numbers of Figs. 1, 2, and 3 for the efficiency 
e vD = (dE vD /dt)/(L„ e +L De +AL Vy ): 



e„j> £ (2...3)10 



(e„ e ) 20 km 



1.5 • 10 52 erg/s 13 MeV R^ sk 



(1) 



where i?disk is the inner radius of the disk. However, 75-80% of 
the annihilation energy is deposited in or near the disk (Fig. 3) 
and will serve to drive a baryonic wind. The fraction of the v 
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energy going into e + e~ pairs in a baryon-free funnel along the 
z-axis and thus being potentially relevant to produce energetic 
7-rays is at most only ~ (0.4-0.9) • 10 50 erg/s at the time of 
maximum luminosity (t <; 6 ms). 

Post-merging evolution. Our simulations suggest that some 
material, possibly ~ 0.1 M , could remain in a disk around 
the central BH. The energy for a 7-ray burst might be provided 
by the v and v emitted from this disk. The efficiency of vv an- 
nihilation increases linearily with the v luminosity (Eq. (1)). 
The v emission properties of the disk should be discussed self- 
consistently with the disk evolution for which the viscous mo- 
mentum transport and viscous heating plays a crucial role. The 
duration of v emission, accretor mass, disk radius, and disk tem- 
perature cannot be chosen and varied independently. We now 
try to relate these quantities by simple conservation arguments 
for disk and v emission. 

The lifetime of the disk will decrease with larger viscosity 
77 because the viscous force that generates a torque carrying 
angular momentum outwards is increased. Setting the rate of 
change of angular momentum equal to the torque exerted by 
the viscous stress yields an estimate for the accretion rate of 
disk material: dM/dt ~ 6irr)h, where 2h ~ 2R & is the ver- 
tical diameter of the thick disk around the BH with mass M 
and Schwarzschild radius jR s ss 9(M/3 M Q )km. This leads to 
an accretion timescale t acc ~ AM i /(6-Kr]R s ) for a disk mass 
AMd. For small 77, the disk will stay cool and nearly transpar- 
ent to v 's. The total v flux L v will therefore increase with the 
viscous energy dissipation rate. The latter is approximated by 
dQ /dt ~ \r\GMjR\ for a (Newtonian) Keplerian disk with an 
inner radius R d ~ 3R S ~ 6GM/c 2 that is taken to be the inner- 
most stable circular orbit around the central accreting BH. The 
timescale of v emission by the disk scales like i a cc oc l/rj and the 
v luminosity like L v oc 77 in the optically thin case. The energy 
converted into e + e~ by vv annihilation is therefore (Eq. (1)) 
E vC , oc L^t dcc oc 77. If the disk viscosity 77 is large, the disk 
will be strongly heated and become very opaque against z/s. 
The v luminosity then behaves like L v ~ Vd(d<5/dt)t acc /i di ff 
for a disk volume Vd, and E v p decreases with 77 according to 

E uD oc l/(r?4ff)- 

The dependence of E vIl on 77 in the ^-transparent and opti- 
cally thick cases suggests that there is a maximum of the annihi- 
lation energy for that value of 77 at which t- dcc « tdiff- Assuming 
the part of the disk where most of the z/'s are emitted to have 
mass AMd and to be a (homogeneous) torus with center at 4i? s 
and radius R s (Mochkovitch et al. 1993, Jaroszyhski 1993), and 
using tdiff ~ 3Rjpa/(m u c) with the v interaction cross section 
a and the atomic mass unit m u , one determines the value of the 
disk viscosity where vv annihilation yields the largest energy 
as 77* ~ |7rm u c/o-. With a ~ 4.3 • 10- 41 (/cT/5 MeV) 2 cm 2 this 
gives 77* ~ 1.6- 10 27 (fcT/5MeV)- 2 gcm- 1 s- 1 . 

Considering a disk with this optimum value 7?*, one can 
estimate the internal temperature T int of the disk from equat- 
ing L v ~ Vd(dQ/df) with L v ~ e v Vi/t m . Here e v ~ 
3 • |a ra dT^ t is the sum of the energy densities of all 3 kinds 
of non-degenerate vv pairs with a ra d being the radiation con- 



stant. One finds (fcT int /5 MeV) « 10(AM d , O i) 1/4 ii Si9 3/4 where 
AM d ,oi = AMj/fO.l M Q ) and i? Sj9 = i? s /(9km) is normal- 
ized to the Schwarzschild radius of a 3 M Q BH. Using this and 
the requirement that L v ~ Si(3^a m aT^ UTl )c/4 (Sd is the surface 
of the toroidal v emitting disk) allows one to determine the 

surface temperature: (fcT surf /5 MeV) « 0.7(AM d ,oi)" 1/8 i?' i / 9 8 . 
With kTim « 50 MeV the optimum value of the disk viscosity 
becomes 77* ~ 1.6 • 10 25 g cm^'s -1 . This corresponds to an 
effective cn-parameter a ss 7)*/(pc s i? s ) ~ (3-6) • 10~ 4 , c s be- 
ing the sound speed and the range of a values accounting for 
differences when the matter is relativistic or non-relativistic. 

For the diffusion and accretion timescales one obtains 
iacc « <diff ~ 0.67(7f)- 1 (AMd,oi) J R s ;9 1 s. This suggests that 
the disk is optically thick for z/s. Under the described con- 
ditions the dynamical shear viscosity due to z/'s is r) v ~ 

0.9 • 10 25 (fcT int /50MeV) 2 (AMd,oi)" 1 -Rs,9g cm_ls_1 which > 
interestingly, is of the right size to account for rf . The v lu- 
minosity is L v ~ 5 • 10 52 (AM d ,oi)" 1/2 ^ i / 9 2 erg/s. The total 
energy radiated away over the time i acc becomes L v t dcc ~ 
3.4 • 10 52 AMd,oi erg which is equal to the gravitational bind- 
ing energy of mass AMd when it is swallowed by the BH 
at r = 3R S : E grdV « GMAM d /(3R s ) = |AM d c 2 . With the 
efficiency of Eq. (1), using L„ e ~ (i ... \)L V ~ 10 52 erg/s, 
(e„ e ) « 3fcT surf ~ 11 MeV, and R^ = 3R S , we find for the 
total energy transformed into e + e~ pairs by annihilation of v 
and v radiated from the disk: 

E vD ~ (2 ... 4) • 10 49 (AM d ,oi) 3/8 R$(* erg . (2) 

This result was derived for the most favorable conditions to 
get large values for E v!l . In particular, it was assumed that the 
dynamical viscosity of the disk attains the optimum value rj*. 
Also, the efficiency e v ^ of Eq. (1) most likely overestimates 
the useful fraction of the annihilation energy by a considerable 
factor. From our hydrodynamical simulations we, moreover, 
deduce that a disk with a mass <; 0. 1 M Q will be formed in NS- 
NS mergers only under very special conditions, e.g., when there 
is high angular momentum in the system due to NS spin like 
in our model C64. The smaller angular momentum in models 
A64 and B64 allows much less material to remain in a disk. 
These arguments might also be true for NS-BH mergers. Finally, 
general relativistic effects were completely neglected in the 
presented considerations. Jaroszyhski (1993) showed that they 
significantly lower the energy that can be transported to infinity. 

Therefore it seems to be extremely hard to account for the 
energies of cosmological 7-ray bursts by vv annihilation. The 
same conclusion was arrived at on grounds of different argu- 
ments by Jaroszyhski (1995) who found that only BH's with a 
high angular momentum (a ~ 1) yield promising models. This 
favors collapsed cores of rapidly rotating WR stars in the "failed 
supernova" scenario (Woosley 1993). Certainly, our analytic 
considerations do not include effects due to possible convective 
overturn or instabilities in the disk which might increase the v 
fluxes and the average v energies considerably and thus might 
help vv annihilation. Also, magnetic fields in the merging NSs 
and in the disk were neglected. Within a lifetime of several 
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tenths of a second initial S-fields of ~ 10 —10 G might be 
amplified by factors of a few 100 in the rapidly rotating disk 
around the BH (rotation periods ~ 1 ms) and might become 
energetically important (Rees, personal communication). 
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